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Low-temperature spin dynamics of the double-layered perovskite La2-2a;Sri+2iMn207 (LSM0327) 
was systematically studied in a wide hole concentration range (0.3 < x < 0.5). The spin- wave 
dispersion, which is almost perfectly two-dimensional, has two branches due to a coupling between 
layers within a double-layer. Each branch exhibits a characteristic intensity oscillation along the 
out-of-plane direction. We found that the in-plane spin stiffness constant and the gap between the 
two branches strongly depend on x. By fitting to calculated dispersion relations and cross sections 
assuming Heisenberg models, we have obtained the in-plane (Jy), iniro-bilayer (J±) and inier-bilayer 
(J') exchange interactions at each x. At x = 0.30, Jy = —4 meV and J± = —5 meV, namely almost 
isotropic and ferromagnetic. Upon increasing x, J± rapidly approaches zero while |Jm| increases 
slightly, indicating an enhancement of the planar magnetic anisotropy. At x — 0.48, Jy reaches 
—9 meV, while J± turns to +1 meV indicating an antiferromagnetic interaction. Such a drastic 
change of the exchange interactions can be ascribed to the change of the relative stability of the 
d x 2_ y 2 and d 3z 2_ r 2 orbital states upon doping. However, a simple linear combination of the two 
states results in an orbital state with an orthorhombic symmetry, which is inconsistent with the 
I<i/mmm tetragonal symmetry of the crystal structure. We thus propose that an "orbital liquid" 
state realizes in LSM0327, where the charge distribution symmetry is kept tetragonal around each 
Mn site. Orbital liquid states are formulated in a theoretical model which takes into account strong 
electron correlation. The calculated results satisfactorily explain the systematic changes of the 
exchange interactions in LSM0327 observed in the experiments. 



I. INTRODUCTION 

Rare-earth doped Mn perovskite oxide Ri-^A^MnOa 
(R: rare-earth ion, A: alkaline-earth ion) is a prototype of 
the colossal magnetoresistance (CMR) materials. Con- 
siderable amounts of efforts have been made to clarify 
the maguetic, electrical and structural properties of these 
systemsEij It has become recognized that the complex 
and delicate balance among the internal degrees of free- 
dom of electrons, i.e., charge, spin and orbital, is a key to 
understand the physics of these materials. M In the Mn 
perovskite oxides, the orbital degree of freedom arises 
from an electron in the doubly degenerated e g states of 
a Mn 3+ ion in a MnOe octahedron. Comparing with the 



charge and spin degrees of freedom, however, the orbital 
degree of freedom has been much less explored, partially 
due to lack of experimental techniques to directly measure 
the ordering processes and the states, with the exception 
of a pioneering work by Akimitsu and Ito who have es- 
tablished the orbital ordered state in K2CUF4 by mea- 
suring the anisotroay of the magnetic form factor using 
polarized neutrons Recently, resonant x-ray scattering 
techniques were successfully applied to the detection of 
the orbital ordering process in Lao.sSri^MnO^ which 
was a significant step toward the understanding of the 
"third" degree of freedom of electrons. In the case of 
LaMnOa, the resonant x-ray scattering techniques have 
provided direct evidence of the orbital ordering.E£!l Two- 



1 



dimensional (2D) planar ferromagnetic coupling found in 
this three-dimensional (3D) lattiea-ean be naturally ex- 
plained by this orbital ordering jlifllJ 

Ri-zA^MnOa consists of MnC>6 octahedra which are 
three dimensionally connected with sharing the corners. 
Due to mismatch of ionic radii of Mn and (Ri_ x A x ) 
ions, the Mn-0-Mn_jpond angle deviates from 180°, which 
is called buckling.E3 Since doping holes inevitably alter 
the average ionic radius of (Ri- X A X ) ion, the amount 
of buckling alsa changes, resulting in a variety of crys- 
tal structures.^ Structural phase transitions are also ob- 
served with changing temperature, a£2 and can be in- 
duced by an external magnetic field. lj The complexity 
of structural properties of Ri_ x A x MnC>3 is an interest- 
ing issue,E2l but makes it difficult to study the role and 
significance of the orbital degree of freedom in affecting 
the magnetic and transport properties, because orbitals 
are strongly affected by surrounding structural environ- 
ment. 

Mn perovskite oxides are generally represented by the 
Ruddlesden-Popper notation (R,A) ra+1 Mn„03„_|_i, where 
the effective dimensionality can be adjusted by changing 
the number of MnC>2 sheets, n, blocked with (R,A)2 02 
layers. As for the, double layered Mn perovskite (n = 2), 
Moritomo et alBS have found an extremely large magne- 
toresistance around Curie temperature T c in a single crys- 
tal of Lai^Sri.sM^Cv, which was followed by intensive 
studies of La9_ 2 xSri +2a; Mn 2 07 (LSM0327) and related 
compounds.EJ Figure |l| schematically shows the strucz. 
ture and the magnetically ordered state at x = 0.40.L2 
A comprehensive magnetic and structural phase diagram 
of LSM0327 has been established by Kubota et a/SEa 
in a wide range of x (0.30 < x < 0.50) as shown in 
Fig. H(a) through systematic powder neutron-diffraction 
studies combined with the Rietveld analysis. They have 
found that LSM0327 exhibits a planar ferromagnetic 
(FM) structure FM-I in the range 0.32 < x < 0.38 at 
low temperatures and that a finite canting angle between 
neighboring layers starts appearing around x ~ 0.39 and 
reaches 180° (AFM-I, i.e., A-type AFM) for x > 0.48. 
They also found that the magnetic moments are aligned 
parallel to the c-axis at x = 0.30, indicating a phase 
boundary between x = 0.30 and 0.32. At x = 0.50, the 
magnetic structure exhibits complicated temperature de- 
pendence due to charge ordering.Ejo 

In contrast to the rich magnetic phase diagram, 
the structure of LSM0327 is fairly simple. Although 
anomalous structural behaviors were reported around T c 
suggesting . a strong coupling among charge, spin and 
lattice,E3t3 there is only a single tetragonal (I4/mmm) 
phase in the entire hole concentration (0.30 < x < 1.0) 
and temperature ranges (T < 400 K) studied so far, ex- 
cept a recently discovered orthorhombic (Immm) phase 
which-exists in a limited concentration range (0.75 < x < 
0.95) Bj This is most likely due to the layered structure 
which absorbs the changes of the Mn-0 bond lengths and 
the average ionic radius of La and Sr ions upon doping. 
It was also found that the Jahn- Teller (JT) type lat- 



tice distortion Ajt of Mn-Og octahedra monotpnically 
changes with increasing x as shown in Fig. 0(c) £3 Note 
that Ajt is defined by the ratio of the averaged apical 
Mn-0 bond length to the equatorial Mn-0 bond length, 
i.e., Aj T = (d Mn -o(i) + ^Mn-o(2))/2dMn-o(3), where 
^Mn-o is a distance between nearest neighbor (NN) Mn 
and O ions. The positions of the O ions are depicted 
in Fig. ||(c). The results are La-good agreement with x- 
ray diffraction measurements. z3c3 The JT distortion in 
LSM0327 stabilizes either d x 2_ y 2 or d 3z 2_ r 2 state. The 
contraction of MnOg octahedron upon doping implies the 
stabilization of d x 2_ y 2 state, i.e., a pseudo-2D e g band, 
in a heavily doped region. As discussed later, such a 
structural change itself cannot account for the systematic 
changes of the exchange interactions within the conven- 
tional double exchange scenario. The dominance of the 
A-type AFM structure with the decrease of Ajt is as- 
cribed to thc_change in the e g orbital state from d 3z 2_ r 2 
to d x 2_ y 2 .E3a The importance of the e g orbital state was 
also-pointed out in striction measurements by Kimup. et 
alEj as well as Argyriou et altB and Medarde et aZ.Ell 

The dynamical magnetic properties of LSM0327 with 
x = 0.4 were measured by Fujioka et alE-i They found 
that the spin wave dispersion is almost perfectly 2D with 
the in-plane spin stiffness constant D ~ 151 meVA. This 
value is similar to that of Lai-^Sr^MnOa (LSM0113) 
with x ~ 0.3, though Tc is three times higher. They 
found that there exist two branches due to a coupling be- 
tween layers within a double-layer. They have analyzed 
the spin-wave dispersion and the differential scattering 
cross section by applying the Holstein-Primakoff trans- 
formation to a Heisenberg Hamiltonian with in-plane 
(J||) and mira-bilayer (JjJ interactions (See Fig. |l|). 
They have estimated that the mfra-bilayer coupling is 
~ 30 % of the in-plane coupling, which is contrary to 
the fact that the Mn-0 bond lengths are similar. They 
speculated that d x 2_ y 2 orbital is dominant at x = 0.40, 
which enhances the double-exchange, i.e., ferromagnetic, 
interaction within the planes. This interpretation is con- 
sistent rWjith the conclusion drawn by previous structural 
studiesE3cj The Z-dependence of the scattering intensity 
due to the spin wave show an excellent agreement with 
the theoretical calculation of the differential scattering 
cross section. Similar inelastic neutron-scattering ex- 
periments were independently performed by Chatterji et 
aZ.,Ej£3 which gives consistent results with that of Fu- 
jioka et aZ.E-3 

The importance of the e g orbital state in determin- 
ing the magnetic and transport properties of LSM0327 
is clear. Moreover, its simple structure makes LSM0327 
more favorable platform to study the roles of orbital de- 
grees of freedom than LSM0113. The resonant x-ray 
scattering technique, however, is not directly applicable 
to LSM0327. This is because the orbital state is pre- 
sumably not antiferro-type long-range ordering as seen 
in LaMn03, where the two types of the orbital are al- 
ternately aligned, thus the superlattice reflections due 
to the orbital ordering do not appear in the resonant 
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x-ray scattering. Instead, it is necessary to determine 
the e g orbital polarization. In the present study, we 
have carried out a series of inelastic neutron-scattering 
measurements on single crystals of the layered perovskite 
La 2 _2xSri +2a; Mn207 (LSM0327) at x = 0.30,0.35,0.40, 
and 0.48. To quantitatively determine the magnetic in- 
teractions in LSM0327, we have calculated the spin- 
wave dispersion and the differential scattering cross 
section numerically by applying the Holstein-Primakoff 
transformation!!!! and the Bogoliubov transformation to 
a Hcisenberg model with the in-plane {J\\), mtra-bilayer 
(Ji) and inier-bilayer (J') interactions. We found that 
the exchange interactions systematically changes with 
changing x. Such an x dependence of the exchange in- 
teractions is well explained by an orbital liquid picture. 



II. EXPERIMENTAL PROCEDURES 

LSM0327 powder was prepared by solid-state reaction 
using prescribed amount of pre-dried La203 (99.9 %), 
Mn 3 4 (99.9 %) and SrC0 3 (99.99 %). The powder mix- 
ture was calcined in the air for 3 days at 1250 -1400°C 
with frequent grinding. The calcined powder was then 
pressed into a rod and heated at 1450°C for 24 h. Single 
crystals were melt-grown in flowing 100 % O2 in a float- 
ing zone optical image furnace with a travelling speed of 
15 mm/h. We powderized a part of single crystals and 
performed x-ray diffraction, which shows no indication of 
impurities. Some of the crystals were also checked with 
electron probe microanalysis (EPMA), which indicates 
no particular spatial inhomogeneity within the instru- 
mental error. The x — 0.30 sample was examined by in- 
ductively coupled plasma (ICP) analysis, which revealed 
that the ratio of La, Sr and Mn is 28.6 : 32.2: 39.2, which 
is in good agreement with the ideal ratio, 28.0: 32.0: 
40.0. These analyses indicate that our samples are suf- 
ficiently stoichiometric and homogeneous. All the single 
crystals studied are domain samples with mosaic spread 
of 0.3 - 0.8° full-width-at-half-maximum (FWHM). The 
samples have a cylindrical shape, which size is typically 
5 4> x 20 — 30 mm. Similarly grown samples in the range 
0.3 < x < 0.5 were powderized and studied in detail 
by powder neutron diffraction techniques and the Ri- 
etveld analysis. The results were already published in 
Refs. HJ24. 



The transport properties of the samples 
were also measured, part of which were published in 
Refs. [38],[39| The results are consistent with previously 
reported dataB3 

Neutron scattering measurements were carried out 
using the Tohoku University triple-axis spectrometer 
TOPAN located in the JRR-3M reactor of the Japan 
Atomic Energy Research Institute (JAERI). The spec- 
trometer was set up in the standard triple-axis mode with 
the fixed final energy at 14.7 meV and the horizontal col- 
limation of Blank-60'-S-60'-Blank. The (0 2) reflection 
of pyrolytic graphite (PG) was used to monochromatc 



and analyze the neutron beam, together with a PG fil- 
ter to eliminate higher order contamination. The sample 
was mounted in an Al can so as to give the (h /) zone 
in the tetragonal IA/mmm notation. The Al can was 
then attached to the cold finger of a closed-cycle He gas 
refrigerator. All the data were taken at 10 K. 



III. RESULTS 

As shown by Fujioka et a/.@ and Chatterji et a/.0@ 
the spin wave dispersion of LSM0327 with x = 0.4 
should have two modes, i.e., acoustic (A) and optical 
(O) branches, due to a coupling between layers within 
a double-layer. It was theoretically shown that the A- 
branch has maximum intensity at / = 5n (n: integer), 
while the phase of the ,-Q-branch is shifted by ir in the 
double layered system Ea We thus measured the spin- 
wave dispersions along [h 0] around (1 0) and (1 5) 
for the A-branch and around (1 2.5) and (1 7.5) for 
the O-branch. To study the differential cross sections of 
spin waves, we have also measured the /-dependence of 
the spin-wave intensities of A and O branches at a fixed 
transfer energy AE = Ei — Ef. 

Figure |3|(a) shows the dispersion relations of spin waves 
at 10 K for x = 0.30. Error bars correspond to the 
full- width at half-maximum (FWHM) of peak profiles in- 
cluding the instrumental resolution. Spin waves of the 
A-branch are well defined in the low q and low energy 
region. However, the O-branch exhibits a large broaden- 
ing even at the magnetic zone center. The /-dependence 
of the constant energy scan at AE = 20 meV are shown 
in Fig. |](a). As expected, the A and O branches exhibit 
intensity maxima at I = 5n and / = |(27i + 1). Solid 
and dotted curves are fitting to theoretical calculations, 
which is described in the next section. 

Figure [|(b) shows the dispersion relations for x = 0.35. 
The dispersion curves of both the A and O branches be- 
come slightly steeper than those of x = 0.30, indicat- 
ing that the in-plane magnetic interaction Jy increases 
only gradually. However, the gap between the A and 
O branches becomes almost half of that at x = 0.30. 
Since the gap corresponds to the out-of-plane magnetic 
interaction Jj_, this result indicates that J± decreases 
considerably. More quantitative analysis will be made 
in the following sections. We have also noticed that spin 
waves are fairly well defined below 20 meV, and that they 
becomes significantly broad above 20 meV, even in the 
constant energy scans. The /-dependence of the constant 
energy scan at AE =11 meV are shown in Fig. |J(b). 

Figure ||(c) shows the dispersion relations foM = 0.40. 
Part of the data has been already reported.^ The l- 
dependence of the constant energy scan at AE = 5 meV 
are shown in Fig. |](c). Following the tendency between 
x = 0.30 and 0.35, the dispersion curves become steeper 
and the gap becomes smaller. As shown in Fig. ||(c), the 
energy- width shows anomalous broadening near the zone 
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boundary, ssihich was found by Fujioka et alxi Furukawa 
and HirotaEll investigated this broadening from both the- 
oretical and experimental point of view, and ascribed it 
to a strong magnon-phonon coupling. Let us consider a 
dispersionless optical phonon branch at HQo, and a spin 
wave dispersion Hu>(q). When a magnon with momentum 
q has energy Lo(q) > Oo, it is possible to find an inelastic 
channel to decay into a magnon-phonon pair with mo- 
mentum q' and q — q 1 ', respectively, which satisfies the 
energy conservation law, u)(q) = oj(q') + Slo • This decay 
channel gives rise to an abrupt broadening of the line- 
width of the spin wave branch which crosses the optical 
phonon. 

Figure ||(d) shows the dispersion relations for x = 0.48. 
The ^-dependence of the constant energy scan at AE — 
5 meV are shown in Fig. [|(d). 

The two branches measured at (1 5) and (1 2.5) 
nearly degenerate. Unlike x — 0.30, 0.35 and 0.40, x = 
0.48 has the A-type AFM (AFM-I) structure, resulting in 
a fundamental difference in the spin wave dispersion. We 
will discuss this difference in detail in the next section. 

As pointed out by Furukawa and Hirota, there exists 
an optical phonon branch around AE = 20 meV for 
x = 0.40, which we have also confirmed for the other 
compositions we studied in the present work. We have 
noticed that the line-width of spin waves above this char- 
acteristic energy of 20 meV become significantly broaden, 
which is consistent with the strong magnon-phonon cop* 
pling model mentioned above. Khaliullin and Kiliancil 
considered an orbitally degenerate double-exchange sys- 
tem coupled to Jahn- Teller active phonons, which ex- 
plains the softening of spin waves at the zone boundary 
found in various ferromagnetic manganese oxides. Their 
model could be applicable to the anomalous broadening 
of the spin waves of LSM0327 near the zone boundary. 
However, the large line-width of the O branch near the 
zone center may not be accounted for because their the- 
oretical model does not affect the small-momentum spin 
dynamics. We thus believe that there exists a signifi- 
cantly strong magnon-phonon coupling as suggested. In 
the present paper, we have combined constant q and en- 
ergy scans to efficiently measure the dispersion relations, 
which is our principle target of the present work. Con- 
stant energy scans are particularly useful to avoid con- 
tamination from dispersionless optical phonon branches. 
To further investigate this issue, however, it is necessary 
to measure the energy widths at various q utilizing con- 
stant q scans, which we plan to carry out in the next 
step. 



IV. ANALYSIS 

In order to analyze the experimental results of the spin- 
wave dispersion relation and the scattering cross section 
in LSM0327, we start from the Heisenberg model where 
3d electrons in a Mn ion are treated as localized spins. 



In the FM (A-type AFM) structure for LSM0327 with 
x =0.3, 0.35 and 0.4 (x =0.48), a magnetic unit cell in- 
cludes two (four) Mn ions termed A and B (A, B, C, 
and D). Between Mn ions, three kinds of the exchange 
interactions, i.e. the in-plane (J\\), intra-bilayer (Ji) and 
inter-bilayer (J') exchange interactions are introduced. 
The schematic picture is shown in Fig. [j]. The Hamilto- 
nian is given by 

+jx & x (?« + + Jl E + (i) 

5± 8' ' 

where S l (fu) is the spin operator at a Mn ion I in the 
i-th unit cell and fu is a position of the ion. Spin quan- 
tum number is assumed to be S = 2(1 — x) + ^x with a 
hole concentration x. Zj_ = I' = (B, A) for I — (A, B) 
in the FM structure and l± = (B,D,A,C) and V = 
(D, C, B, A) for / = (A, B, C, D) in the A-type AFM one. 

8\\ , S± and 8' indicate the vectors connecting the nearest 

neighboring (NN) Mn ions; ^ = (±a,0,0) and (0, ±a,0) 
where a is a distance between NN Mn ions in the ab plane. 
8 X = (0, 0, c x ) and 8' = (±o/2, ±o/2, -d) for I = B and 
D, and 8± = (0,0, -cx) and 8' = (±a/2,±a/2,cf) for 
I = A and C, where c± and d are distances between 
NN MnC>2 layers and NN bilayers, respectively. By ap- 
plying the Holstein-Primakoff transformation to Eq. (Q), 
the Hamiltonian is rewritten as follows, 



(2) 



Here, ip(k) = (ak,bk) for the FM structure and ip(k) = 
(ofc, b\, Cfc, d\) for the A-type AFM structure, at, bk, Cfe 
and dk are the boson operators for the spin operators S A , 
S B , S c and S D , respectively. s(k) is given by 



e(k) 



x y 
y* x 



(3) 



with 



x = — 4J||5|l — i(cosa/s x + cosafc y )| 

~J ± S - AJ'S, (4) 

and 

y = ,J ± Se lk * c± +4J'5cos 



(5) 



for the FM structure, and 



s(k) = 




(6) 
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with 



X = — 4J||5|l — -(cos ak x + cos ak y ) 
+JiS - AJ'S, 



ik z c± 



and 



Y = J ± Se 



Z = 4j'Scos{^)co S (^)e^ 



(7) 



(8) 



(9) 



for the A-type AFM structure. By utilizing the canonical 
transformation of the Hamiltonian Eq. (^), the dispersion 
relations of the spin waves are obtained. As for the FM 
structure, in particular, the dispersion relations are given 
analytically aswt = x± |y|. On an equal footing, the dif- 
ferential scattering cross section for the inelastic-neutron 
scattering from spin wave is given by 



dQdcu' ~ mc 2 \2 9 W> 



2 k> 



k 8 



X E E ^ - 5 @ m + K)U} m (q)U mV (q) 

U'm qG ^ 



+ 6(u> + <)<S(Q -G + qjn^U^ml^q) \, 



(10) 



where q and uj™ are the momentum and energy of spin 
wave of the mode m, respectively, and n™ = l/{eP u i —1) 
is a Bose factor with temperature T = 1/(3. U(q) is a 
matrix introduced in the canonical transformation, Q = 
ki — kf is the momentum transfer with ki (kf) being the 

momentum of the incident (scattered) neutron, and G 
is a reciprocal lattice vector. F(Q) and W(Q) are the 
magnetic structure factor and the Debye- Waller factor, 
respectively. 

The experimental results of the dispersion relation and 
the differential scattering cross section in LSM0327 are 
fitted by utilizing the least-squares method. The calcu- 
lated results are shown in Figs. || and || together with 
the experimental data. Note that the dispersion rela- 
tion along [h 0] and the cross section along [0 I] are 
sensitive to Jii and J±, and J± and J', respectively. In 
x = 0.3, 0.35 and 0.4, the A and O branches are well sep- 
arated. These two correspond to the in-phase and out-of 
phase motions of spins in NN MnC>2 layers. The energy 
separation between the two branches at the point T and 
the stiffness constant of the A-branch in the ab plane are 
given by — 2S(J± + 4 J') and D = —SJ\\, respectively. 
An intensity oscillation along [0 I] is factorized by the 
functions l + cos(c_Lg z ) and 1 — cos(cj_g z ) for the A and O 
branches, respectively, where c± ~ c/5 with c being the 
lattice constant in the c axis. This is attributed to the 
spin correlation between NN MnC>2 layers controlled by 
J± . Additional fine structures in the intensity are caused 



by J' . On the contrary, in the A-type AFM structure, 
four modes of the spin wave exist and separate into the A 
and O branches corresponding to the in-phase and out-of 
phase motions of spins between NN bilayers. Here, the 
each branch is doubly degenerate and the energy sep- 
aration between the two is of the order of J'. Since 
J' ~ J||/1000 at x = 0.48, as mentioned later, these 
cannot be observed separately by the experiments. An 
intensity oscillation along [0 1] is factorized by a func- 
tion 1 — cos(c^g 2 ) originating from the antifcrromagnetic 
spin alignment between the NN Mn02 layers. 

The x dependence of the exchange interactions is 
shown in Fig. ||. All interactions systematically change 
with x; with increasing x from 0.3, |Jm| increases, Jj_ 
rapidly approaches to zero and changes its sign from neg- 
ative to positive. | J'\ decreases with x, although its value 
is 10^1000 times smaller than \J» | and | J±\. The system- 
atic change of the interactions correlates with that of the 
lattice distortion in a MnOg octahedron represented by 
Ajt (Fig. ||(c)). However, this is an opposite direction 
predicted by the conventional double-exchange scenario, 
where |J||| is reduced with increasing x because magni- 
tude of the double-exchange interaction is proportional 
to the hopping integral between Mn ions in the strong 
Hund-coupling limit. Therefore, a structural change it- 
self cannot account for that of the exchange interactions. 
Let us taking into account the orbital degree of freedom 
in a Mn ion. In LSM0327 with hole concentration x, 
1 — x electrons occupy the two e g orbitals. Character of 
the occupied orbital controls the anisotropy of the hop- 
ping integral of electrons, i.e., that of the ferromagnetic 
double exchange interaction. The systematic change of 
the exchange interactions can be explained by assuming 
that the id x 2_ y i orbital is relatively stabilized with in- 
creasing x. The exchange interaction between Mn ions 
is sum of the ferromagnetic double-exchange interaction 
and the antiferromagnetic superexchange one J afm act- 
ing between ti g spins. The more the d x 2_ y 2 orbital is 
stabilized, the more the ferromagnetic interaction in ab 
plane (along the c axis) becomes strong (weak). Then 
J afm overcomes the ferromagnetic interaction along the 
c axis as shown in the region of x = 0.4 — 0.48. This 
assumption for th e orbit al stability is consistent with the 
x dependent A,/Tcirc2lO and is supported by the previ- 
ous theoretical work where the stability of the orbitak is 
examined by the Madelung potential calculation.OE^ 



V. SUMMARY AND DISCUSSION 

We discuss possible orbital states in LSM0327 with 
0.3 < x < 0.5 and its relation to the anisotropy of the fer- 
romagnetic interaction in more detail. The orbital state 
at each Mn ion is represented by the pseudo-spin opera- 
tor defined by 



(11) 
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for /i = (x,z). di-y S is the annihilation operator for the 
e g electron at site i with spin s and orbital 7, and are 
the Pauli matrices. In the eigen state of T iz = +(— )l/2, 
an electron occupies the d 3z 2_ r 2 (d x 2_ y 2) orbital at site i. 
Tiz (Tix) describes the charge quadrupolc moment with 
tetragonal (orthorhombic) symmetries and couples with 
the lattice distortion with the same symmetry; 



(12) 



lfl — X,Z 



where Qi Z (Qi X ) describes the O ion distortions in a 
MnOg octahedron. The orbital ordered state is char- 
acterized by a magnitude and an angle of this operator, 
i.e. \(f)\ = V^) 2 + <T Z ) 2 and 9 = tarr^T,)/^)) 
where (• • •} is the thermal average. For example, in the 
(d^ x 2_ r 2^ G? 3 j,2_ r 2)-type orbital ordered state observed in 
LaMnC>3, O = 27r/3 and — 2tt/3 for the Mn sites where 
the d^ x 2_ r 2 and d^ y 2_ r 2 orbitals are occupied, respec- 
tively. 

As mentioned in the previous section, a relative weight 
of the occupied d x 2_ y 2 orbital increases continuously with 
increasing x from 0.3 to 0.48 in LSM0327 where the crys- 
tal structure remains to be tetragonal (I4/mmm). That 
is, (T z (k — 0)} is gradually reduced with keeping the 

condition (T x (k = 0)> = where f(k) = l/NJ2i e i% ' ?i Ti 
with the number of the Mn ion N and the position of 
the i-th ion fi. This cannot be satisfied by the uniform 
orbital ordered state where one kind of orbital character- 
ized by O is occupied at all Mn sites. This is because the 
change of the orbital state is represented by the rotation 
of (T) in the (T Z )-(T X ) plane. One may think that the 
antifcrro-type orbital ordered state explains the experi- 
ments, when the condition = — 0^, with Qa(b) being 
the angle in the orbital space for the A (B) sublattice, 
is satisfied. However, this is ruled out by the experi- 
mental fact that the expected superlattice reflection was 
not reported ip -0.3 < x < 0.5 by the x-ray and electron 
diffractions MScJ 

One of the possible orbital states realized in LSM0327 
is an orbital liquid state. This state was originally pro- 
posed in the ferromagnetic metallic state in LSM0113 



by Ref. 46, where the orderings of both T z and T x are 



suppressed by the low dimensional character of the or- 
bital fluctuation. In the case of LSM0327, (T z ) is finite 
due to the layered crystal structure. On the other hand, 
Ti x does not show ordering and symmetry of the charge 
distribution remains to be tetragonal at each Mn site. In 
order to formulate this orbital state, let us start from the 
Hamiltonian where the intra-site Coulomb interactions 
in Mn ions are taken into account; 

-J„J2Si-S ti + Aj2Tiz, (13) 

i i 

where di ia = di ia (l — rii l9 ){l — n,;^ (J )(l — n^) is the an- 



nihilation operator of an e g electron excluding the doubly 
occupied states of electrons due to the strong Coulomb 
interaction. Si is the spin operator for an e g electron 
defined by Si = | J] ss ' 7 d\ ls a ss 'di 7S < and S ti is the spin 
operator for ti g electrons with S = 3/2. The first and 
second terms in Eq. (|l3|) represent the hopping of e g elec- 
trons between NN Mn sites and the Hund coupling be- 
tween e g and t2 g spins, respectively. In the third term, 
A describes the splitting of the energy levels of d 3z 2_ r 2 
and d x 2_ y 2 orbitals due to the tetragonal distortion Ajt 
of a MnO@ octahedron. It is shown from the theoreti- 
cal calculation in Ref. |4^ that A monotonically decreases 
with increasing x for LSM0327 and its maximum value 
is of the order of 0.5 eV. Instead of the actual crystal 
structure of LSM0327, a pair of the 2D sheets where a 
squared lattice consists of the Mn ions is introduced be- 
cause of the weak inter-bilayer exchange interaction. We 
adopt the slave-boson scheme where di ia is decomposed 
into a product of operators: di 1<y = //t^s^ where fi 
and Sia- are bosonic operators for charge and spin de- 
grees of freedom, respectively, and Ti 7 is a fermionic one 



for orbital associated with the constraint of V„ s]„ s 

L — < a i(T 



Y,y T n T H and fifi + Ea 4r S «7 = 1 at each site - The 

mean field approximation is introduced; (fjfj) = x and 

E ff ( s lv) = { l ~ x ) e v withe^ = +(-)l for a ferromag- 
netic (antiferromagnetic) bond in the orbital part of the 
mean-field Hamiltonian. It is well known that the slave- 
boson mean-field approximation is suitable to describe 
the spin liquid stateMU The ratio of the ferromagnetic ex- 
change interaction in the ab plane to that in the c axis is 
given by R = J\\/J± = Xab/Xc with xi = T,j-y'( T l 7 T i+h') 
for I — ab and c. The results are shown in Fig. |^(a) 
where Jh is assumed to be infinite and io is the hopping 
integral between d 3z 2_ r 2 orbitals in the c axis. R contin- 
uously decreases with increasing A implying that d x 2_ y 2 
orbital becomes stable relatively. This feature does not 
depend on the hole concentration in the calculation. 
Since A is expected to continuously decreases with x in 
LSM0327, the results in Fig. 0(a) explain the experimen- 
tal results presented in Fig. g(b). A dependence of {T x ) 
are shown in Fig. ^(b). We note that (T z ) is zero. It 
is clearly shown that the continuous change of R is con- 
trolled by the character of the occupied orbital. We also 
present the schematic pictures of the spatial distribution 
of the electronic charge at a Mn site for the proposed 
orbital liquid state in Fig. [?]. The charge distributions 
have tetragonal symmetry and are not represented by 
any linear combination of the atomic wave functions of 
the d 3z 2_ r 2 and d x 2_ y 2 orbitals. 

To summarize, we have systematically studied low- 
temperature spin dynamics of the double-layered per- 
ovskite La 2 _ 2a; Sri +22; Mn207 (0.3 < x < 0.5). The 
acoustic and optical branch of the 2D spin-wave disper- 
sion relations as well as characteristic intensity oscilla- 
tions along the out-of-plane direction are successfully ex- 
plained by theoretical calculations assuming the Heisen- 



G 



berg models with the in-plane (J\\), mira-bilayer (J±) 
and mier-bilayer (J') exchange interactions. We have 
found that the ratio R = J\\/J± drastically decreases 
upon doping holes, which indicates that the d x 2_ y 2 or- 
bital becomes more stable than the d 3z 2_ r 2 orbital. Since 
a simple linear combination of the two states results in 
an orbital state with an orthorhombic symmetry, incon- 
sistent with the IA/mmm tetragonal symmetry, we have 
introduced an "orbital liquid" state, in which the charge 
distribution symmetry is kept tetragonal around each Mn 
site. 
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FIG. 2. (a) Structural and— magnetic phase diagram of 
LSM0327 after Kubota et ozH (b) Several different mag- 
netic structures appearing in the phase diagram are schemat- 
ically drawn, (c) Hole concentration dependence of the JT 
distortion, which is defined as the ratio of the averaged apical 
and the equatorial Mn-0 bond lengths, at room temperature 
and 10 K. 



FIG. 1. Schematic representation of the magnetic spin ar- 
rangement on Mn ions in the IA/mmm tetragonal cell of 
Lai.2Sri.gMn207. Each Mn ions are surrounded by an 06 
octahedron. The lattice parameters are a = b — 3.87 and 
c = 20.1 A at 10 K.B. Notations are explained in §. |fv|. 
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FIG. 3. The dispersion relations of spin waves at 10 K for 
(a) x = 0.30, (b) x = 0.35, (c) x = 0.40, and (d) x = 0.48. 
Error bars correspond to the FWHM of peak profiles. Solid 
circles and open circles indicate the acoustic branch and the 
optical branch, respectively. Solid and dotted curves are ob- 
tained by fitting to theoretical models described in §. pv| . 



FIG. 4. Differential cross sections along the out-of-plane 
direction, which were obtained from the Z-dependence of the 
constant-i5 scans. The solid and open circles indicate inten- 
sities of the acoustic and optical branches, respectively. Solid 
and dotted curves are obtained by fitting to theoretical mod- 
els described in §. [iv|. The acoustic branch is dominant at 
(10 0) and (1 5), and the optical branch is dominant at 
(1 2.5) and (1 7.5). 
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FIG. 5. (a) x dependence of the exchange interactions ob- 
tained by analyses of the dispersion relations and the scat- 
tering cross sections in the inelastic neutron scattering ex- 
periments, (b) x dependence of the ratio of the exchange 
interactions J± and J« . 
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FIG. 6. (a) Theoretical results of the exchange interac- 
tions ratio R=J±/J||. J\\ and J± are the exchange interac- 
tions between NN Mn sites in the ab plane and along the c 
direction, respectively, (b) The relative number of occupied 
electrons in two e g orbitals 2{T Z ). 




FIG. 7. Schematic pictures of the charge distribution for 
e g electrons in the case of (a) A = —0.4, (b) and (c) 0.4. 
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